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SUMMARY 


A general Fokker-Planck equation is deduced, which describes the 
distribution of geomagnetic ally trapped electrons as a function of longitude, 
time, energy and mirror point field intensity. A special variable for the 
longitudinal position of a particle is introduced. 

The physical interpretation of this equation is analyzed for several 
special cases. In particular, it is found that the usual procedure of aver- 
aging over longitude in order to obtain a longitude -independent equation, 
is not valid for the description of electrons mirroring at low altitudes. In 
that case, the longitude dependence cannot be averaged out, and a four- 
dimensional equation must be used. 

The coefficients representing longitudinal drift, ionization loss and 
multiple, screened Coulomb scattering in the general Fokker-Planck equa- 
tion are deduced. An approximation is given for a simple model of the 
field and the atmosphere in the South American Anomaly. Qualitative 
conclusions are drawn for the mirror point "flow" along lines of force, 
for a stationary electron distribution drifting through the Anomaly. It is 
concluded that the region East of the Anomaly, initially depleted by pre- 
cipitation, is replenished by electrons whose mirror points were situated 
in a narrow "window" in B, before passing through the Anomaly. A very 
limited extension of the atmosphere in the region of the Anomaly is re- 
sponsible for this mechanism of replenishment. A diurnal effect for the 
electron fluxes in the replenished region is predicted, controlled by the 
diurnal variation of the atmosphere in the Anomaly. 
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calculations of longitude dependence of 

GEOMAGNETICALLY TRAPPED ELECTRON FLUXES 

PART I 

THE GENERAL FOKKER-PLANCK EQUATION 
FOR ELECTRON DIFFUSION 


INTRODUCTION 

Time has come to look with more detail into the latitude dependence of geo- 
magnetically trapped particle fluxes from the theoretical point of view. Until 
recently, experimental information on this subject was very scarce. Only after 
the discovery of the so-called South American Anomaly by Vernov, et al . 1 and 
the more detailed surveys with Discoverer satellites by Mann, Bloom and West 2 , 
attention got more and more focussed on the experimental analysis of particle 
precipitation in the Anomaly, and the mechanism of subsequent replenishment 
of the depleted shell regions. In particular, Imhof and Smith 3 , Paulikas and 
Freden 4 and Mihalov et al. s have made a careful study of artifically injected 
electron fluxes, measured by various satellites. On the other hand, Freden and 
Paulikas 6 and White 7 have analyzed proton fluxes at low altitudes in the Anomaly. 
Evidence for Bremsstrahlung X-rays from electrons precipitating into the 
Anomaly was recently found by Ghielmetti, et. al. 8 . It is, therefore, desirable 
to set up a theoretical description of the longitudinal behavior of trapped particles, 
and to test by comparison with experimental data, the various assumptions made 
about interaction processes governing particle diffusion. 

Theoretical description of electron trapping, diffusion and precipitation was 
so far done only for configurations averaged over all longitudes 9 ’ 10,11 ’ 12,13 . 

In these papers, a time dependent Fokker-Planck equation was set up for the 
electron distribution function and used to follow the evolution in time of a given, 
initial electron flux. 

In order to study the longitude dependence of the electron distribution on a 
given magnetic shell, it is necessary to derive a more general diffusion equation 
which contains an additional variable related to longitude. In this equation the 
coefficients, which depend upon the atmosphere, will be a function of longitude 
and local time of day in the Anomaly. One important by-product of the equation 
derivation will be the identification of an appropriate variable by which to describe 
the longitudinal dependence. We will call this variable X. It will be used, together 
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with the particle energy E , the scalar magnetic field B at the mirror point, and 
the well-known Mcllwain shell parameter L, 14 to describe the space in which 
particle densities change with time,t. 

The main purpose of Part I is to set up such a general equation, to discuss 
its physical meaning, to compare it with the previously used longitude-averaged 
equation, and to draw some general, qualitative conclusions about longitude de- 
pendence of trapped electrons. Part II* will deal with atmosphere-field configura- 
tions to be used in this longitude dependent description. In Part III*, results of 
a numerical integration of the general equation will be presented. 

Before setting up our equation, let us picture the problem in general terms. 
Consider the familiar B-L space, in which trapped radiation fluxes are usually 
described (Fig. 1). It can be shown that electrons whose mirror points are below 
100 km cannot remain trapped for more than a few bounces. Thus, it is qualita- 
tively useful to consider the 100 km level as the location of a sink. We shall 
return to this point in more detail later. At a given longitude the 100 km level 
can be displayed as a locus B c (L,X) in B-L space. Essentially no electrons 
will be found with mirror points above this curve, that is, for mirror point fields 
larger than B c (L, X) . We have plotted two extreme loci of the 100 km level, 
corresponding to a longitude (that is, a value of X) right in the "center" of the 
South American Anomaly, and a longitude over the Pacific Ocean, respectively. 
From Fig. 1, one can see that for particles of a given energy, the 100 km curve, 
B c (L, X), "oscillates" with longitudinal drift frequency between these two extreme 
positions. As B c (L, X) lowers from its maximum, or "Pacific", position, the 
region in B-L space between the two extreme positions, called the "shadow 
region", is wiped clean of particles. As B c (L, X) rises from its minimum, or 
"Anomaly", position, the opportunity exists for atmospheric scattering, or any 
other nonadiabatic process, to repopulate this shadow region. This wiping clean 
and repopulation has been likened to the action of a windshield wiper. 

Satellite observations do, in fact, find particles in the shadow region (Imhof 
and Smith 3 ). This observation means that particles have enough non-adiabatic 
interaction during one drift around the earth to diffuse into the wiped-out or 
"shadow" region. Our main goal is to study this process of fast replenishment 
of the shadow region in the case of electrons interacting with the atmosphere by 
Coulomb scattering and energy loss. 


Will be published at a later date. 




Figure 1-The “windshield wiper" effect in B-L space. 


REPRESENTATION OF THE ELECTRON DISTRIBUTION 
AS A FUNCTION OF LONGITUDE 

We shall describe the trapped electrons by the distribution function used by 
Welch, Kaufmann and Hess 10 . For a given L-shell we denote by 


8N = U(B, E, «p, t) 80 SB SE 


( 1 ) 


the number of electrons contained at the time t in a tube of field lines of mag- 
netic flux 80 situated at a longitude (for the time being, we shall use the geo- 
graphic longitude of the magnetic equatorial point), with mirror points between 
B and B + 8B, and with kinetic energies between E and E + SE. We shall neglect 
L-scattering; this is why we do not include explicitly L as a variable in U. 

The distribution U is not directly measurable. It is, however, the physical 
magnitude which most appropriately describes trapped particles. Relations with 
experimentally accessible quantities are given in 10 . In particular, the number 
Sn of particles contained in the volume element 8ASs of a tube at a point where 
the field is B' , and which mirror in the interval dB at B > B', (Fig. 2) is given by 
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B 



Figure 2 


Sn = U(B, E, <p, t) dB SE B SA 


6s (2) 

V| ( (B', B, E) r b (B, E, cp) 

The ratio in (2) represents the probability of 
finding a particle of the population (1) in the 
volume element SASs. r b (B, E, cp) is the half- 
period of bouncing of these particles, i.e., the 
time it takes to go from one mirror point to 
its conjugate. In the real geomagnetic field, 
it depends slightly on longitude. v| | (B\ B, E) 

= v(E) i/l - B'/B is the particle's velocity 
parallel to the field line, at the point B' . For 
a dipole field, the total path of a particle dur- 
ing a bounce period is vr b ^ 2LR (R = radius 
of earth). The total particle density per unit 
energy at a given point B' is therefore 


n(B' , E, cp, t) 


Sn 

SA Ss SE 


b* r B c 

V (E) L 


U(B, E, <p, t) dB 


r (B.E, cp) i/l - B'/B 2RL J, /I _B'/B 


b' r ° 

2RL 

*'B 


UdB 


(3) 


The integration is performed along the field line between the B value at the given 
point, and some upper cut-off B c , beyond which all particles are absorbed in the 
atmosphere. The density (3) is related to the omnidirectional counting rate C of 
a spherical detector by 


C = 



vn A dE 


(4) 


where A(E) is the effective cross sectional area of the counter for particles of 
energy E . 

Let us now consider electrons of a given energy, mirroring at the same B , 
trapped between two neighboring shells labeled I and I + SI (Fig. 3). I is the 
second adiabatic invariant 14 of these particles computed for field lines belonging 
to the inner shell: 
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I + 81 is the value of the second invariant, taken along a line of force belonging 
to the outer shell, between mirror points with the same field intensity B. B 0 is 
the field at the equatorial point of a line of force (point of minimum B on a given 
shell); cp is the geographic longitude of this equatorial point. We call % I = 
n Sl/Sy the gradient of I at the equatorial point of a magnetic shell (Figs. 3,4). 
Northrop and Tellgr 15 have shown that the equatorial drift velocity, averaged 
over one bounce, is given by 


u o = 


m v -> — 

V I x B = -2- -2- b 

B^ t. 0 eB„ 


eK 6 T b 


'0 T b 


(5) 


where p is the particle's momentum, e the charge (<0 for electrons) and r b is 
the half-period of bouncing. For the different vectors, see Fig. 4. 



Figure 3— Shell geometry and elementary flux tube 
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Figure 4-Unit Vectors, mean equatorial drift 
velocity and gradient of I 


We denote with Sx = u Q St the distance between equatorial points of two 
neighboring field lines of a given shell; 8x is then the element of arc of the 
shell's equatorial B 0 -ring; x =/dx is the total arc length along the equator from 
a given initial point and can be used as a label to locate field lines on a given shell. 
Evidently 


and 



Sep 


( 6 ) 



(7) 


where dx/3<p is a field-geometric factor (Bx/Bqp = 27r/360 v RLfor a centered dipole), 
and ep the geographic longitude of the equatorial point. Let us take a tube of lines 
of force of equatorial cross section Sx 8y and magnetic flux B 0 §x8y (Fig. 3). By 
the definition of U (1), the number of electrons in this flux tube, which mirror 
between B and B + SB and with kinetic energy between E and E + SE , is given by 


SN = UB 0 Sx Sy SBSE 


( 8 ) 
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We now want to follow the history of these particles, as a function of time. 
Before setting up the Boltzmann-type equation governing the distribution U, we 
have to introduce a convenient variable for the longitude, i.e. a convenient label 
for the line of force of a given magnetic shell, around which a particle is instan- 
taneously spiralling. We cannot take 9 orx (7) as suitable variables, because 
U is a distribution in flux , not in 9 or x . And in the flux expression B 0 Sx 5y 
intervening in (8), not only Sx, but also Sy is a function of longitude, for the 
general geomagnetic field. 

In order to find the correct longitudinal variable, we have to transform the 
flux into an expression Sep = B Q SXSY, in which SY - still related to the radial 
increment Sy — is now longitude-independent, so that it may be ignored when 
following the particles during their longitudinal drift. In that case, X would be 
the correct longitudinal variable, its differential SX containing complete infor- 
mation about the longitude dependence of the flux S 4> of a tube filled with particles, 
as they drift around the earth. 

In order to find X and Y, we just have to determine the longitude dependence 
of Sy, the equatorial distance between two neighboring magnetic shells. Notice 
that this dependence is a purely geometric one, independent of the particles’ 
dynamic variables. Let us take Mcllwain's definition of L, 14 which we write 
here: 


r 3 l 3 b f /i 3 b\ 

M \ M j 


( 9 ) 


(R: radius of the Earth, L dimensionless, I in cm). We further consider the 
relation: 


B 


0 


M 

R 3 L 3 


( 10 ) 


(only valid within about a few percent, because L actually fluctuates along a 
line of force 14 - 16 ). Combining (9) and (10), and differentiating, we have 



3I 2 


F' SI 


R 3 L 3 


( 11 ) 
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for a constant mirror field B. F' is the derivative of Mcllwain's function with 
respect to its argument. Dividing by Sy, we obtain, for absolute values 


si_ _ B o § y r 2 l 2 

Sy 3/RL j2 p» 


( 12 ) 


Notice that 1/B 0 SB Q /Sy 1/B 0 V 0 B as Sy 0, where V Q B is the gradient of 
|B| taken at the equatorial point. For a current-free magnetic field, 1/B 0 V 1 B 
is equal to the curvature of the line of force at the point where the normal 
gradient is taken; in our case, 1/B 0 V Q B is therefore the curvature at the equa- 
torial point. We shall introduce the dimensionless number 


K 


0 



3/RL 


(13) 


which is the curvature of a line of force at the equatorial point, in units of the 
equatorial curvature 3/RL of the corresponding dipole line (k 0 = 1 for all dipole 
shells). In the general case, k q depends on longitude. From (12) and (13), we 
finally have: 


where 


SI 

Sy 


= 7.1 =— 

0 /(B) 


/(B) = 


I 2 F' 


r 2 l 2 


(14) 


(15) 


is a dimensionless function of the mirror point field intensity B. Notice that, 
according to (5) and (14), we have 


u 


o 


P K o 
eB o r b^ 


(16) 
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and for the angular drift velocity, taking into account (6): 
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1 

3x U ° 


3<p 


1 p K o 

Bx eB o A. 

Bq) 


( 17 ) 


It is important to point out that A.(B) (15) is a "pure dipole function", even in the 
case of a slightly distorted dipole as the real geomagnetic field (this is not true 
for T b ). It is easy to verify that our function X is identical with 1/6E, where E 
is the function calculated by Hamlin et al. 17 (formula 22 in this reference). All 
this of course is valid only within a few percent (fluctuation of L on a magnetic 
shell). Notice finally that k 0 contributes to the longitude dependence of the drift 
velocity. However, one must not forget that r b also depends on longitude in the 
real field. According to (14), 


Sy = 


KSI 


(18) 


In this relation, the only longitude dependent quantity is k 0 . We then introduce 
the variables Y and X in the following way: 


SY = k 0 Sy = \SI a Y = \I 


(19a) 




(19b) 


X increases from West to East. Notice that labeling field lines by the coordinate 
X is equivalent to picking up lines of a shell with a "density" (number per unit 
equatorial ring arc length) inversely proportional to the curvature k q (13). For 
a pure dipole field, or for the outer shells in the real field, « 0 = 1, and X is 
identical with x (7). For low-L shells, k 0 may differ from 1 as much as ±6% 

(See Part II). With the new coordinates (19a) and (19b), the flux S 4> is 


Sqb = B 0 Sx Sy = B 0 SX SY 
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in which SY is now independent of longitude. Adopting X (19b) as the appropriate 
longitudinal variable, we can write the number of electrons in a flux tube as 


SN = U(B, E, X, t) SB SE SX • B Q SY 


( 20 ) 


As B 0 SY is a constant throughout the whole forthcoming discussion we shall drop 
it at once whenever we consider the number of particles in a flux tube. 

Finally, it is important to remark that there is a more exact expression for 
the fundamental relation (14), which takes into account the fact that in the real 
geomagnetic field, L differs from 1/R (M/B 0 ) 1/3 (10) by a small, longitude 
dependent amount AL = AL(L, B, cp ) 16 . This relation is given by 



1 

X(B) 



However, the correction term RL 3/3y(AL/L) is expected to be small as compared 
with k 0 , and the longitude variations of k 0 . On the other hand, the correction term 
implies the existence of a shell splitting as a function of the mirror point field 
intensity 16 . But in the whole treatment which follows, this splitting is being 
neglected, and all electrons, initially on a given line of force, are supposed to 
populate the same shell, regardless of their mirror points. 


DERIVATION OF THE FOKKER-PLANCK EQUATION 

FOR LONGITUDE DEPENDENCE 

We are now in condition to set up the general equation governing the distri- 
bution of trapped electrons. These electrons will undergo displacements in B, 

E and longitude, caused by three types of mutually independent interactions: 

(1) A change in the mirror point field B , due to the stochastic process of 
multiple Coulomb scattering in the atmosphere; 

(2) A change in kinetic energy E , due to ionization slowing -down in the 
atmosphere, considered here as a non -stochastic process (i.e. neglecting 
straggling); 

(3) A change in longitude due to interaction with the static magnetic field 
(longitudinal drift), again a non- stochastic process. 
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All these processes are physically independent of each other, although all inter- 
vening parameters are in general functions of the three variables B, E and 
longitude. 

Let us relate the distribution of electrons as it appears in (20), with the 
distribution of electrons at a slightly earlier time t - At, at a different longi- 
tudinal position. For the time being, we shall forget the non-stochastic character 
of E and X, and treat all variables as if they were of the same, stochastic nature. 
We write: 

U(B, E, X, t) = 


# <B -A E-e, X-A t-At) nfa-A E-e, X-<f, A e, A At) d/3ded£ +Q(B, E, X, t)At 

( 21 ) 

In this relation, the distribution of electrons which at the time t are at a position 
X , is linked to the distribution of those electrons which at an earlier time t - 
At were at X ~A and which happened to scatter, slow down and drift the right 
amount in the interval At, in order to become part of the population described 
by the left hand of (21). 

ndySde d£ is the a priori probability that these electrons have undergone 
just the right changes in mirror point field, energy and longitudinal position, in 
the time interval At. 

The source term QAt represents the contribution of electrons added by 
injection to the original bunch of particles, during At . 

Expanding all intervening functions in Taylor series in /? , e , g and At , we 
obtain: 


# (U <»> = - i (U < £ » - s A s <u </?>) + 


+ Q + higher order terms in £ 2 , e 2 , e/3, etc. 


( 22 ) 
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The brackets stand for the average per unit time of the enclosed variables: 



A II d/3 de d£ 


We can now re-instate the non-stochastic character to E and X , by taking 11 
as a delta function in energy and longitude: 

ri(B, E, X, A e, A At) = n 0 (B, E, X, A At) S (e + k At) S (£ - <f At) 

In this expression, e, is the change of energy per unit time, due to ionization 
loss. Evidently, 


e= — {e} (23) 

T b 

where {e} is the energy loss in one half-bounce (path from one mirror point to 
its conjugate). We have taken {e} and e as positive quantities, and evidenced 
the energy decrease by the sign in the first delta function. 

As to the drift velocity in the X coordinate, we have, according to (19b), 
(16) and (17): 


i = E, X) 


dX u o _ 1 3x . p 1 

dt Kq K 0 Bcp 9 e B 0 T b \ 


(24) 


n 0 is the probability for a change in the mirror point field due to Coulomb 
scattering, in the interval At. It is evidently given by 


n 0 (B, E, X, A At) = P(B, E, X, P) - (25) 

r b (B, t!', A) 

where P is the probability for B -scattering during one half-bounce. 

In all this we have implicitly assumed that scattering and slowing down are 
extremely small during one half-bounce of the electron. In other words, we 
suppose that P << 1 for fi f 0 and approximately = 1 for /3 = 0 (this imposes 


12 



a limiting condition to the applicability of our discussions when dealing with 
very high atmospheric densities). 

Equation (22) finally becomes, up to the second order in /3: 


ft + = i *7 ^ (u » + Q 


(26) 


where 


{/3> {/3}=J/3Pd/3 

' 7 T b 

wi th 

<f/3 2 )> = — {/ 3 2 } i/3 2 } = J /3 2 P d/3 (27) 

W T b 


The angular brackets are average changes per unit time; the curled brackets 
{ } represent average changes per half-bounce. 

We can re-write the longitude-convection term as a function of ordinary 
geographic longitude, taking into account (24): 



— — ( U — — cp 

1 3x Bcp l k q Bcp 

K 0 3< iP 


Bx eB Q X. 


Bcp 



(28) 


Notice that in a pure, centered dipole, * = 1, Bx/B<p and cp constants, and the 

convection term would reduce to cp BU/Bcp. 


Equation (26) is the most general Fokker-Planck equation describing the 
longitude and time dependence of a trapped particle flux, which scatters and 
slows down in the atmosphere. 


The initial condition for equation (26) is the distribution U given at a certain 
longitude, for a given initial time. There is a "natural" boundary condition at 
the magnetic equator: 


^ = 0 for B = B 0 (29) 
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At the ends of lines of force in the high atmospheric density region, there is 
no physical boundary condition for U as a function of B : as particles enter higher 
densities, they get slowed down by ionization and die away in energy space . 
Therefore, it would be incorrect to impose a condition of the type U = 0 for 
B t B c . Moreover, equation (26) is not valid at all near such fictious boundary 
(for derivatives and integrals could not be interchanged any more when passing 
from (21) to (22)). Of course, if one is only interested in electron fluxes in B 
regions which always remain high up, i.e. far away from high atmospheric den- 
sities, a boundary condition for U as a function of B may well be adopted for 
practical, computational reasons. 


DISCUSSION OF THE GENERAL EQUATION 

Let us now discuss equation (26) from the point of view of its physical mean- 
ing. First of all, in absence of scattering, slowing down and sources, the right 
hand of equation (26) is zero. We then have: 


~+ (U £) = 0 (30) 

dt dX 

This equation tells us that a given initial distribution of electrons injected at 
a time t 0 at a line labeled X 0 , will proceed drifting eastwards with speed t , 
always keeping the electron distribution inversely proportional to the local drift 
speed £ : 


u£ 


cons t 


(31) 


As t is energy-dependent, equation (30) represents a "flight-time spectrometer 
effect" on the particles, after injection. Notice carefully that expression (31) 
means constancy along the drift path of a given bunch of particles. It is not a 
constancy in time. 


If we now integrate equation (26) over one complete longitudinal cycle of X, 
for a fixed time t , the second term on the left vanishes. We obtain: 



Ue 




dX + 0 


(32) 
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v( c 


Dividing byj> dX and calling 


U 


AV 


1 UdX 

f d4s 



J *0 


r 


(33a) 

<t> dX 

T — 


J 

J K o 


i eUdX 




£ AV =' 


UdX 


(33b) 


f</ 


AV 


n 


,3/UdX 


| UdX 

| <^>U dX 




AV 


I 


UdX 


(33c) 


we obtain the longitude -independent equation 


BU. 


AV B 


a, - 3E <°» 4 «> - 5 b i u » v 9 1T l u *» W.V I * 5 (34 > 


1 3 2 * 


This equation is formally equivalent to the time -dependent equation used by 
several authors 9 ’ 10 ’ 11 - 13 . Notice, however, the following remarks: 

(1) The average distribution U AV is computed integrating over the new 
variable X, which according to (33a), means that U must be weighted 
at the different longitudes with the inverse of the equatorial curvature 
of the field lines. 

(2) The "coefficients" e AV , and y are not simple averages over 

the new coordinate, but they are weighted with the distribution U itself. 

In other words, they are functionals of the unknown distribution. 
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All this leads us to the conclusion that a longitude-average Fokker -Planck 
treatment of the problem of trapped particle diffusion is not valid , unless we 
restrict our description to only those particles which mirror high enough at all 
longitudes so that their distribution function can be expected a priori as very 
little longitude dependent. Equation (34) is certainly meaningless as a simple 
differential equation for that portion of B-L space, which descends below about 
300 - 350 Km in the South American Anomaly. This is precisely the domain 
where the longitude -independent treatment has so far failed to give numerical 
results compatible with experimental measurements. 

In order to "legalize" the longitude -independent description for such cases 
in which we know a priori that U will not depend strongly on X , we have to evaluate 
the coefficients (33) for a distribution U which has the smallest possible longi- 
tude dependence . This dependence is given precisely by (31), i.e. for the case 
of absence of interactions with the atmosphere. Introducing (31) in (33), and 
taking into account (24), we obtain "true" coefficients (independent of U): 



These expressions are formally identical with the longitude-average coefficients 
used by Hassitt 18 , in which the atmospheric constituents (the only strongly longi- 
tude dependent variables actually contained in e , (^fi ) and ) are weighted 
inversely proportional to the equatorial drift velocity at each longitude. 

We must, however, insist again that the use of a longitude -independent 
Fokker-Planck equation for the description of trapped electron diffusion, is not 
valid at all for those high B values, for which the atmospheric interactions in the 
region of the Anomaly will introduce a notable departure of U from the "collision- 
less" expression (31). 
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Let us now return to equation (26), for the steady state case BU/3t - 0. 
This is very nearly true for natural radiation belt electrons during quiet solar 
wind conditions, and even for artifically injected fluxes, provided we are inter- 
ested in the longitude dependence of the electron distribution, at times after in- 
jection long compared to a typical longitudinal drift period. In that case, we 
can also neglect the contribution of the natural source term Q . Taking into 
account (24), (23) and (27), we introduce the distribution function 


w = = u -E- _L 

eB o T tA 


(36) 


and the coefficients 




{e} 



— 1 / \ e B. 

£2 = ± (r2\ = k _JL {/? 2 } 

£ p 


which now represent changes per unit longitudinal coordinate X . Notice that 
the factor 





P 



= b(B, E) 


(38) 


physically represents the "bounce-density," i.e., the number of half-bounces per 
unit longitudinal coordinate, which a particle of mirror point field B and energy 
E makes at any position on the shell. Notice that this bounce density is indepen- 
dent of the longitudinal position X . 
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Using (36) and (37), the steady state form of (26) reads: 


3W 

dX 



2 dB 2 


(W/3 2 ) +^r( We) 

ok 


( 39 ) 


This is the fundamental equation which describes the longitude dependence of a 
stationary, source free electron distribution, now represented by W (36). It 
will be the main subject for the rest of this work. 


A useful concept for a qualitative analysis of equation (39) is that of the 
"mirror point flow in B-X space," i.e., the mean change of the mirror point 
field intensity of a particle, per unit longitudinal coordinate. Let us consider 
the following relation ~ 


/•B(X) 

I WdB = const, (40) 


where the integral is taken along a given field line between the equatorial point 
B 0 and a generic field point B(X) . We shall call the relation B = B(X) , which 
ensures the constancy of (40), "characteristic trajectory of a mirror point" (see 
Fig. 5). Notice carefully that, in general, B = B(X) is not the real path of the 



Figure 5 
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mirror point of a given electron, in B-X space. However, it does become the 
actual average path, if we neglect energy loss and also neglect the longitude 
denendence of r. . In that case. is nrnnnrtionni to the nnmher of nartieles 

* L> ‘ \ - L- - L - - - - - k~ 

mirroring between equator and B(X). But even if these conditions are not satis- 
fied, the function B = B(X) defined in (40) is a useful concept. 

In absence of all interactions, B(X) = const., i.e. the mirror point field 
intensities remain unchanged. We now turn on Coulomb scattering, but still 
neglect energy loss. We then introduce the quantity v m = dB(X)/dX as the 
"mirror point flow in B-X space," obtaining, with the use of (39) and (40): 


_ dB(X) I f B(X) il dB = -= 1 * BlnW 

dX ~ W J B BX P ' 2 3B ‘ 2 P BB 


(41) 


As defined in (41), v m describes the average "flow" along the field lines, of 
mirror points of electrons of a given energy, as they drift in longitude. If we 
wish to add the energy loss mechanism, an additional term would appear in the 
right side of (41), and particle individuality would be lost. 


Notice that there are three physically quite distinct contributions to the 
mirror point "flow." The first term of the right hand of (41) represents a steady 
increase of B , (ft is positive), i.e. a steady lowering of mirror point altitude. 
This contribution eomes_from the first order, "streaming" term in B (39). In 
absence of dispersion (ft 2 =0), /3 is the only contribution to v m ; in that case, 

B = B(X) = Jv m dX is just the characteristic of the corresponding first order 
differential equation. 


The second and third terms on the right hand side of (41) arise in the dis- 
persion mechanism. Their presence in (41) clearly shows that dispersion also 
contributes to a steady, average "flow" of mirror points along the lines of force. 
These two terms are governed by the gradients (along field lines) of the coeffi- 
cient ft 2 , and of W (or U), respectively. The gradient of ft 2 is always positive, 
so that the second term in (41) always represents an upward motion of the par- 
ticles's mirror points (back-scattering from the denser atmosphere). If, on the 
other hand BW/3B < 0 , the third term represents a lowering of the mirror point 
altitude. If W / 3 2 = const, along a line of force, the contribution from the two 
dispersion terms is zero: The number of mirror points thrown upwards by the 
gradient in ft 2 , exactly compensates the number of mirror points streaming 


v'atvi n i v» o 


downwards due to the gradient of W . The only mirror point flow which iviuaiiia, 
in this case, is that of the first order term: v = ft. 
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DETERMINATION OF THE COEFFICIENTS 

In this section we shall obtain expressions for the various coefficients inter- 
vening in (39). According to (37), the principal quantities to be found are k , {e}, 
{/3} and {/3 2 }. 

With respect to the dimensionaless function k , defined in (15), we recall 
that it is inversely proportional to a dipole-type function computed by Hamlin, 
et. al. 17 A good approximation of k for the interval 1/4 < B 0 /B £ 1 is given by 


\(B) 2 0.384 + 0.128 exp 


1 B 
1.37 B 0 


(42) 


The energy loss per half-bounce 


{ e } is given by the integral 


{e} 


ps(B) 


3E 

3s 


ds 

» B 


(43) 


The origin of the field line arc length is taken at the geomagnetic equator. 
ds//l - B’( s)/B is the element of trajectory of the electron at a point s where 
the field is B' . Using tables given by Berger and Seltzer, 19 we obtain the follow- 
ing very good approximation: 


BE 

3s 


= 10' 20 N ff ion h(E) (kev/cm) 




where 


2.61 + 349.0 E* 0 - 844 for E < 470 kev 


r < 44 > 


h(E) =- 


4.54 


for E ^ 470 kev J 
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and 


N 


f f . : 


N + 0.881 N . + 0.273 N. 

oxygen nitrogen helium 


(45) 


is the "effective" atmospheric number density for ionization loss of electrons. 

As for the time being we are interested in low L shells only, we did not include 
the contribution from atmospheric ions and electrons. 

In order to evaluate {/ 3 } » the average change of B per half-bounce, we first 
have to perform the average of collisons over the isotropic azimuthal distribution 
of scattering angles. Following Welch, Kaufmann and Hess 10 and slightly changing 
their notation, we obtain 


? = -j)si n2 » (46) 

/3 is the change in mirror point field intensity B , when the electron's pitch 
angle scatters an amount 6 at a field position B' , averaged over all possible 
azimuthal angles of scattering. Now we have to find the average of sin 2 0 at the 
given field point B'. This will be given by 

( 47 > 


sin 2 9 66 (48) 

d9 

is the average contribution from one atom of class i . N. are the number den- 
sities of the different atmospheric constituents at the field point B' ; do- /dd are 
the differential cross sections for screened Coulomb scattering, for each con- 
stituent. Using Moliere's expression for dcr./d<9, 20 and following closely Welch, 
Kaufmann, and Hess, 10 we obtain 


where 


sin 2 9 


) *ir/2 
0 
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where 


N 


e f f . s ca 1 1 . 


= N + 0.774 N + 0.058 N. .. 

oxygen nitrogen helium 


(50) 


is the "effective" atmospheric number density for screened Coulomb scattering 
of electrons. Again, the contribution of atmospheric ions and electrons is 
excluded. 

We finally have to perform the integration along the path of a particle from 
one mirror point to its conjugate: 



Again following Welch, Kaufman and Hess, 10 we obtain 






sin 2 <9 


(52) 


Therefore: 


{/S 2 } = 2B 2 




(53) 
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If we now call: 


S X (B, X) = 


/•s(B) 

*'-s(B) 


B ^ef f scatt (s, X) j 

B'(s) VI - B'(s)/B 


(cm -2 ) 


S 2 (B,X) 


ps(B) 

J-sCB) 


N eff scatt ( S> ^ 

Vl - B'(s)/B 


(cm -2 ) 


SUB, X) 


1 


S (B) 


ef f 


(s.X) 


.(B) /I - B'(s)/B 


ds (cm -2 ) 


(54) 


K(B, E) = 4 x 10 -22 B 



(gauss cm 2 ) 


C(E) = 10 -20 h(E) 


(kev cm 2 ) 


we can summarize expressions (51), (53) and (43) in the form 
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(55) 


(py =kL-Is\ 

{/ 3 2 } = 2BK (Sj - S 2 ) 

{e} = CS' 


DISCUSSION OF ELECTRON DRIFT AND SCATTERING 
THROUGH A MODEL FOR THE ANOMALY 

Before setting up numerical integrations for low altitude electrons, it is 
very helpful to discuss Equation (39) qualitatively, using a simple, approximate 
model of the atmosphere along field lines, in the region of the South American 
Anomaly, and where most of the scattering and energy loss occurs. 

We shall restrict ourselves to the region between 250 and 750 Km altitude, 
and to longitudes within about ±30° of the "center" of the Anomaly (supposed at 
a longitude cp a = 318°E). In the present discussion we shall neglect the energy 
loss term and suppose Bx/Bcp = const, and k 0 = 1. These latter suppositions 
lead to a proportionality between X and the longitude 9 of the equatorial point 
of a field line (19b). 

Using the Harris -Priester atmosphere 21 for low solar activity (S = 70), 
and Southern Hemisphere geomagnetic field B-L rings as calculated by 
Stassinopoulos, 22 we obtain a quite good fit for the effective number density 
(50) as a function of B and 9 , of the type: 


N 


ef f sea 


tt = N o ex P 


GJ - B 

a 

1 AB 


exp 



(56) 


Due to the supposed proportionality between 9 and X, it is easy to convert this 
expression into a function of X . All this is valid only for very low L values. 


24 



For L - 1.25, we have:* 


13 hs local time at the center of the Anomaly: 

N q = 10 11 cm -3 B g = 0.243 gauss AB = 5.25 x 10~ 3 gauss Acp = 4.69° 
02 hs local time at the center: 


N 0 = 10 11 cm" 3 B a = 0.243 gauss AB = 4.23 x 10" 3 gauss A<p = 3.78° 

For an expression of the form (56), with a "scale height" AB very small 
compared to B , the main contribution to the integrals (54) comes from a region 
close to that mirror point which lies in the high atmospheric density end of the 
line of force. We can therefore approximate: 


/B - B \ 

s i = s 2 = N o ex P ) ex P 




3s 


(57) 


and 



(58) 


3B/3s is the average gradient of B along the field line, in the region of interest. 
Taking into account (55), the following relation between the coefficients {/ 3 } and 
{/ 3 2 } holds: 


{/ 3 2 } = 2AB { /3 } 


(59) 


*For higher L shells (L > 1.8), the approximation (56) becomes very crude. Qualitatively, as L 
increases, <p Q shifts to the East, and Acp increases considerably. 
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This is similar to a relation found by Walt and McDonald 12 for the coefficients 
of a time dependent (but longitude -independent) Fokker-Planck equation set up 
in terms of mirror point altitudes, in the region of high atmospheric densities. 
Relation (59) is expected to hold quite generally in the region of high densities, 
provided one inserts for AB the value AB = N , ,/ |3N £r /3B I , where the deriva- 
tion is taken along a line of force. 

With (59) and the linear relationship between X and cp , Equation (39) becomes 


3W _ J_ 3W 

3X ~ 3X Bcp 


~ (W/3) + AB — (W/3) 
H 3B 2 


(60) 


in which fi is of the form 


0 = K ' exp {^r) exp (- 


K' is obtained from (37), (54), (55) and (57), and is a slowly varying function of 
B (slow in comparison with exp (B/AB ). 

The mirror point flow (41) is now 


v 

m 


/3 - AB M _ AB B 
3B H 


91n W 
BB 


Taking into account (61) we can approximate: 


A 

BB “ AB 


(62) 


Therefore 

v ^ - AB B 31n W 
m P 3B 


(63) 


26 



Physically, the cancellation of the first term with the second term in (41) means 
that the upstream of mirror points due to diffusion from the denser atmosphere 
C’backscattering" at lower altitudes) exactly compensates the downstream due 
to the first order term. What is left, is a mirror point streaming, entirely due 
to diffusion, and which may be upwards or downwards, according to the gradient 
of the actual electron distribution, along a line of force. 



Figure 6-Qualitative description of the evolution 
ofa mirror point distribution at very low altitudes. 


In particular, if 3U/3B = 0 for some value B j; there will be no net flow of 
mirror points across this B value. An initial distribution of electrons at a 
longitude cp Q like the one shown in Figure 6, will therefore keep the maximum 
at the same B x , but will broaden towards both sides as the particles drift towards 
the east ( <p 2 > cp x > cp 0 ). This smearing-out will be faster at the high-B side 
(towards lower altitudes), for is much greater there (see expression (63)). 

A similar result was obtained numerically by MacDonald and Walt. 12 

If now W is of the form 


W ~ 


1 


K' exp 



taA \ 
V' 71 / 
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we have a solution of Equation (60), which represents a steady state in longitude. 
In this case, the mirror point flow (63) is approximately equal to Ji : 


v y - AB p 

m ' 


31nW ^ -g 

m ~ p 


(65) 


It can be shown that such a steady state form (64) is attained very soon, provided 
the gradient of W along a line of force is negative. It is also easy to verify that 
(64) is an "equilibrium" configuration: Any small departure of W from (64) will 
cause mirror point flows which will tend to correct such a departure. 

For all these reasons, it is extremely instructive to inspect in detail the 
"motion" of mirror points for such a longitude independent equilibrium state 

(64) . Although, of course, all this is only a crude approximation, we expect all 
qualitative, physical features to be valid also in the real case. This qualitative 
analysis will then help us to set up and interpret the tedious numerical calcula- 
tions (Part HI). 

We shall analyze the characteristic trajectories of mirror points B = B(X) 
(40) throughout the South American Anomaly, described approximately by (56), 
for an equilibrium state electron distribution (64). First of all, we will do our 
discussion in terms of the more familiar longitude variable cp , rather than X, 
which in our description anyway is supposed proportional to the former. In 
order to obtain the characteristic trajectories in B-cp space, we have to integrate 

(65) , taking into account (61) and the relation (19b), which now reads 
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360 


RLcp . 


Starting from an initial longitude cp 0 , well West of the center of the Anomaly, 
and for a shell L = 1.25, we obtain characteristic mirror point trajectories for 
300 kev electrons, schematically represented in Figure 7. Remember that they 
correspond to a distribution of the type (64). In this figure, lines of constant 
altitude for the approximate expression (56) are also shown. Below 250 Km, the 
N eff = N eff (B, cp) dependence is still nearly exponential in B and |<p - <p a | , but 
with a smaller "scale height" AB. We can safely assume for this discussion, 
that particles getting below this level will be removed by precipitation. Above 
750 Km, the scattering effect is extremely small during one drift through the 
region of interest. 
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298 

Figure 7-Characteristic mirror point trajectories in the region of the South American 
Anomaly (model (56)), for a steady state distribution of 300 Kev electrons of the form 
(64), and for L = 1 .25 (13.00 LT in the Anomaly). Electrons whose mirror points are 
initially (i.e. West of the Anomaly), between B . and B , populate the “shadow” 
region. Those between B m<jx and B WE precipitate East of cp a ; those with B > B EW 
precipitate on the West side. Constant altitude lines for h = 250 and 350 Km are shown. 
Energy loss is neglected. 




Notice in Figure 7 that there is a narrow "window" B . - B , at the 

initial longitude cp 0 West of the Anomaly, which contains the mirror points of 
those electrons, which populate the high-B or "shadow" region, East of the 
Anomaly. After one more turn around the earth, all these electrons are precip- 
itated into the dense layers. A certain portion of these (B max - B WE ) will manage 
to get "around" the center of the Anomaly, and precipitate into its Eastern side. 
Electrons mirroring initially at higher altitudes than those corresponding to the 
"window," will lower slightly their mirror points during each passage through 
the Anomaly until they fall into the critical window. Please remember carefully 
that all this corresponds to an equilibrium distribution of the form (64), neglecting 
energy loss, and remember that the mirror point "trajectories" are average 
paths, not actual patterns for a given electron. 

The most striking result of this analysis, is that if multiple Coulomb scatter- 
ing is the only process responsible for the replenishment of the "shadow" region 
East of the Anomaly, then there is a very limited region of the atmosphere in 
the South American Anomaly (limited in both, longitude and altitude), which en- 
tirely governs this replenishing process. We therefore should expect that in the 
real case, electron fluxes in the "shadow" are strongly influenced by the instan- 
taneous atmospheric structure in this critical region. 

In particular, we may expect a diurnal variation of electron fluxes and elec- 
tron precipitation, determined by the diurnal variation of the atmosphere within 
about ±5° of the center of the Anomaly, between about 250-400 Km altitude. We 
conclude, that it would be unrealistic, to carry out numerical calculations with 
a 24-hour averaged atmosphere. Rather, numerical integration should be per- 
formed for different local times at the Anomaly (Part II). 

Notice finally, that according to this picture, replenishment occurs right at, 
or shortly "after," i.e. East, of the Anomaly. Diffusion into the "shadow" region 
far East of the Anomaly is negligible, for Coulomb scattering. We may conclude 
our discussion by pointing out that in the picture described above, there is a 
steady flow of particles out of the lower B region of a shell, towards higher B 
values, and from there, into the atmosphere at the Anomaly. In an equilibrium 
state, this loss must be exactly balanced by a continuous injection from a source 
mainly effective in the lower B region. In absence of such a source, i.e. in a 
non-equilibrium state, this flow and subsequent loss should determine the lifetime 
of the trapped radiation in question. 

As a numerical example, we have calculated the position of the critical 
"window" in B (for the definition of B . and B see Figure 7), for local times 
in the Anomaly as 1300 LT and 0200 LT. The values obtained for L = 1.25 are 
shown in Table 1. Although the B values differ only very little, remember that 
they are in a region where the distribution U is a very steep function of B . 
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Table 1 


1300 LT 
B . 

min 


B 

ma x 

0200 LT 
B . 

min 

B 

max 


300 kev 
0.198 gauss 
0.201 


625 kev 
0.205 gauss 
0.209 


0.207 0.213 

0.210 0.216 


Detailed quantitative results on electron fluxes can be obtained only by num- 
erical integration of Equation (39). In particular, the energy loss term, neglected 
in the above discussion, has to be taken into account; this is particularly impor- 
tant if we analyze the behavior of electrons which dip below about 300 Km altitude 
in the Anomaly. 
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